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Abstract 

The linear response approach to nuclear transport has been extended to pair 
correlations. The latter are treated within a mean field approximation to a pairing 
interaction with constant matrix elements G. The constraint of particle number 
conservation has been accounted for on a time dependent average, which leads to 
modified response functions, both in the pairing degree of freedom as well as in the 
shape variable. The former is expressed by the gap parameter A and the latter by 
a Q which specifies the elongation of a fissioning nucleus. The tensors for friction 
and inertia corresponding to these two collective coordinates are computed along 
the fission path of 224 Th for temperatures around T = 1 MeV and less. It is seen 
that dissipation decreases with decreasing temperature and increasing pairing gap 
and falls well below the values of common " macroscopic models" . Both friction and 
inertia show a sensible dependence on the configurations of the mean field caused 
both by shell effects as well as by avoided crossings of single-particle levels. 
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1 Introduction 



Pair correlations are vital for understanding many elementary features of nuclear physics 
at zero or small thermal excitations. In particular, it is known that they have great 
impact on collective properties, see e.g. [p] or 0. For the case of zero temperature this 
has been demonstrated in various ways. Here it may suffice to mention a few examples 
related to nuclear fission. Since the early 70'ties it is known that pairing will modify 
greatly the effective inertia ||, which in turn may change the penetrability through the 
barrier by orders of magnitude. In addition, the pair degree of freedom in itself may 
portray important features, like pair vibrations etc. [J]. A convenient way of handling such 
properties is by introducing the gap parameter as an additional, independent collective 
degree of freedom, which first has been utilized in J|. Later this parameter has been 
introduced also to the generator coordinate method ||, where the influence of pairing 
vibrations on the spontaneous fission probability has been studied. 

It is of great interest to account for pair correlations also in the description of typical 
transport problems of dissipative systems. On general grounds it is to be expected that 
pairing will greatly diminish nuclear dissipation. Indeed, there exist early formulations 
within linear response theory where such effects have been studied, see [0-0. The present 
approach is meant as a generalization of these in a four fold sense: (i) We will account for 
particle number conservation on a time dependent average by introducing and exploiting 
modified response functions, (ii) We will calculate response tensors (and thus tensors for 
transport coefficients) for a pairing mode plus a shape degree of freedom, (iii) This will 
be done along a complete fission path parameterized in terms of Cassini ovaloids, thus 
improving the results of || for pair correlations, (iv) We will also address temperatures 
below 1 MeV, for which pair correlations become even more important. 

The paper is organized as follows: In Sec. ^| we will quote basic elements of linear 
response theory. For a more detailed exposition of the theory we refer to the review 
article || . Sec.^| contains the derivation of the mean-field Hamiltonian and the response 
functions with pairing included. The Strutinsky renormalization will be applied to the 
calculation of the free energy in Sec. |3.2.1| . In Sec. |3.3| the evaluation of the collisional 
width for quasi-particles will be explained . In Sect.[| the modification of the response 
function will be given which accounts for the average particle number conservation in 
the dynamical process. Finally, in Sec.^ numerical results for response functions and 
transport coefficients are shown, the latter being compared in Sec. |5.3| with those of other 
approaches, with special emphasis put on the temperature dependence of dissipation. 
Section ^ contains a short summary. 



2 Linear response theory for collective motion 

In this section we give a brief review of a description of collective motion within the locally 
harmonic approximation (LHA). It is based on the common hypothesis (see e.g. and 



101 ) of the existence of a set of collective c-number variables which in parametric 
way portray average dynamics of the nucleus as a whole. Fluctuations in these quantities 
shall be neglected; how the latter may be included can be found in 0. The coupling 
of these collective variables to the nucleonic degrees of freedom is traced back to the 
deformed shell model potential which shows up in the (many body) Hamiltonian H sm for 
independent particle motion, to be generalized later to the quasi-particle picture when 



pairing correlation are to be included. 

As is well known, the total energy of the system may be obtained from H sm after 
applying Strutinsky's renormalization procedure. For later purpose it will be convenient 
to account for this fact by adding to the H sm a c-number term £ pot to get 



pot ; 



which in some sense may be called the Hamiltonian of a "renormalized mean field". At 
zero intrinsic excitations, the additional term can be seen to be given by the negative 
value of the average potential energy, see |jj and c.f. HlOj . Such a construction warrants 



that the total energy of the system can be expressed by the mean value (i7 rm /(£;,p;, Q^))- 
It may easily be extended to finite temperature and the correction term may simply be 
found by requiring that the energy obtained from the Strutinsky procedure (at finite T) 
be identical to the one calculated from H rm f, see ||. 

As mentioned earlier, we will be concerned with situations of finite thermal excitations. 
In this case it must be expected that the pure independent particle model breaks down even 
close to the Fermi surface, in the sense that the particles (or quasi-particles) experience 
"collisions". Formally, the latter may be represented by adding to the H rm f a residual 
two-body interaction which leads to the Hamiltonian 

H(xi,pi, Qy) = H rmf (xi,pi, Qp) + V$ s {x h pi) (2) 

However, it is too cumbersome to work with such a form in genuine sense. For this reason 
we are going to approximate the effects of the by dressing the single particle energies 
with self-energies having both real and imaginary parts. A description of details of this 
procedure will be deferred to a later section. We may note here, that this approximation 
shall be done in such a way that these self-energies are insensitive to changes in the 
collective coordinates Q M , for which reason the have been left out in the arguments of 
f( 2 ) 

res 

This latter approximation goes along with the conjecture that the "generators" of 
collective motion are given by the following one-body operators 

a / a . OHixi^puQ^) dH rmf (x h p h Q,,) 

^,P„0,.)- (3) 

These generators make up the coupling between collective and intrinsic degrees of freedom. 
Within the LHA this may be seen as follows, (c.f. Suppose the actual Q M is close to 
some Q ^. One may then expand the effective Hamiltonian to second order as 

H{QM = H{Ql) + - Ql)F, + - - Ql) (Q v (t) - Ql) ( ^— ^ ) 

ii z nv \ u Hpy K 4v I qo to 

"(4) 

to describe general features in the neighborhood of the Q ^. Evidently, such a procedure is 
possible for all static quantities. In the dynamic case one needs to require that the Q^(t) 
does not move away from Q° within some time 6t which is larger than the typical time 
for the nucleonic degrees of freedom, say the one which describes their relaxation to local 
equilibrium. 

The generator F M appearing in (||) is the one of (0) but calculated at It represents 
the nucleonic part of the coupling between nucleonic and collective degrees of freedom 



which is linear in the latter. As implied by @, within the LHA this is the only coupling 
term left. In the second order term the nucleonic part appears only as a static average 
of the corresponding operator. This average is to be built with the density operator 
Pqs(Q°) which in the quasi-static picture is defined by the Hamiltonian at Q Q , namely 
H{Q Q ^). The p qs (Q°) is meant to represent thermal equilibrium at Q ^ with the excitation 
being parameterized by temperature or by entropy. The simplest possibility is offered by 
using the canonical ensemble, or more generally, the grand canonical distribution. For any 
(static) the latter is defined as 

p v {Q^T } fj) = —ex^{-H'(x h Pi,Q^^)/T), Z' = tr(exp(-H'(x i ,p i ,Q ti ,fi)/T)) (5) 
with 

H'(xi,pi, Qp, fj) = H{xi,Pi, Qfj) - fiN (6) 

and /i, N being the chemical potential and the particle number operator, respectively. 

For the remaining part of this section we are going to restrict ourselves to a fixed N 
and discard complications from a possible variation of particle number. Later we are going 
to study the general case for which \i is considered just another "collective variable". 

Before we address time dependent forces let us quote a general relation for the gener- 
alized static ones along some given direction Q^. With E(Q ft , S) being the internal energy 
at fixed entropy S and the T{Q^,T) the free energy at given temperature this relation 
reads 

/dH\ qS _ ( dF(Q„T) \ f dE(Q„S) \ 
W/ { dQ, ) TQ ^ \ dQ v ) S>Q ^ 

It is valid at any both at or away from global minimum where these derivatives vanish 
TJJ and H, again). 



sec 



The equations of average motion for Q^(t) may be constructed looking at energy 
conservation. For an arbitrary path through the multi-dimensional collective space the 
change of the total energy can be expressed as 

It must vanish if we picture the nucleus as a whole as an isolated system. Requiring this 
condition to hold true for any path through the collective landscape one gets the set of 
equations 

To put them into convenient form one first needs to express the average as a 

functional of Q^t) as well as a function of possible changes in either temperature or 
entropy. The form of the coupling term as given in (Q) invites one to apply linear response 
theory to get the actual time dependent part. Considering Yl^F^Q^t) — QV) as a time 
dependent perturbation one may express the deviation of the (F fl ) t from their quasi-static 
values as (see [|]) 

W = -E / x^-s){QAs)-Ql)ds (10) 

,, J — oo 



where x^ v is the causal response function 



X fl At-s) = Q(t-s)Ui (p qs (Q ,T )[F'(t),Fl(s)})^2 l e(t-s)x''(t-s) (11) 



Here, the time dependence of F^(t) (in interaction representation, if one wishes) is deter- 
mined by the H{Q Q ^}. The Fourier transform of ( PH| ) reads 

= - E X^Q u {u) with 8Q v {u) = Q v {u) - Q^M (12) 

V 

with 5Q v {uS) being the Fourier transform of (Q„(s) — 

Piecing together all relevant parts, in Fourier representation and exploiting vector 
notation the set of equations of motion may finally be written in the following compact 
form 

5(F) u = k- 1 t l (u) with q = Q-Q m (13) 

and where the Q™ defines the center of the local oscillators (see |§ where all details have 
been worked out carefully for the one-dimensional case). In ( |T3"D a tensor fc M „ appears 
whose elements often are interpreted as coupling constants, see below. At zero excitation 
it would simply be given by the negative value of the average of the second derivatives of 
the Hamiltonian flQ| . At finite excitations additional terms show up which involve static 



susceptibilities. Which ones will appear depends on the situation one chooses. In || it 
has been argued in favor of working at constant entropy. It is much simpler, though, to 
assume that temperature does not change within the time lapse 5t for which the equations 
of collective motion are derived. For the sake of simplicity we will adhere to the latter 
case, knowing that in this respect we are along the lines of most of the papers in this field. 
For such a situation the coupling tensor, which often is expressed by its inverse k = k ~ , 
is given by (c.f. ||) 



K 



^-\^rwr) +x^(o)-x£,= ^ m +x^(o) (14) 



Here, i s ^ ne static response and 

d(F,r 



T _ ^\ r f 



dQ u 



(15) 



the isothermal susceptibility defined as the change of the quasi-static expectation value 
(F fl ) qs with Q v at given T. It may be noted that at lower temperature the difference 
;w(0) — xliv — known to vanish exactly at T = — may often be neglected as compared 
to the first term on the right of (|I4D. Indeed, such a case will be encountered below. 

Finally, we aim at equations of motion for the of a structure similar to that of 
Newton's equation with a dissipative force. To be able identifying uniquely all forces, as 
well as the corresponding coefficients for inertia, dissipation and local stiffness, it is useful 
to introduce a coupling U ext (t) to some "external" fields f^ xt {t) which finally will show 
up as in-homogeneous terms. A convenient form is given by 

u ext (t) = Y J hr;\t) (i6) 



Defining a response tensor X co u(u;) by 

*W, = -ECH/fH (it) 



it can be shown (see JTTJ and |J, for instance, or |IJ for the undamped case at T = 0) to 
attain the form 

XcoiiM = «(k + xM) _1 xH (18) 

Its pole structure is determined by the collective excitations of the system, whence the 
name "collective" response tensor comes from. 

To finally be able introducing the transport coefficients for average collective motion, 
as expressed in terms of the q^t), one first rewrites ( fP7D in the form 

["Xcoii ~ V)«]<?M = -«r» = -<T» (19) 

Here, (|13j) has been used and the external fields q ext for the g^-modes are defined through 
the last equation on the very right. The set ([19]) turns into the conventional form of 
Newtonian dynamics 

J2(M^q„(t) + 1 , v q v {t) + C^it)) = -qf(t) (20) 

V 

if one may approximate the quantity kx co h _1 (w)k by a second order polynomial in fre- 
quency 

( K Xcoii -M^J 2 - i-f^u + (21) 

Evidently, the coefficients M M „, 7^ and stand for the elements of the tensors for mass, 
friction and stiffness. 

Without any doubt, the transition ( pl|) is crucial to the question as to which extent 
Markovian type of equations may be justified for the very complex situation of nuclear 
collective motion. In the past, two version have been worked out for this transition (see || 
for a review and a guide to original work). On the one hand, one may study the collective 
strength distribution as given by the dissipative part of the collective response. Whenever 
this distribution has prominent peaks (for the mult i- dimensional case see |H]|) they may be 



fitted by Lorentzians, which in turn may be interpreted as the oscillator response function 
associated to the form (|20|) of the classical equations of motion. Another possibility is 
offered by simply expanding the functions given on the left hand side of ([H]) to second 
order in uj around uj = 0. In the one-dimensional case one gets, see ||), 

X(0) + C(0) 



c 



7 



i ^xcoiiM)- 1 ( x (o) + c(o)) 2 



k 2 duo 



uj=0 



X 2 (0) 



C(0) (22) 
7(0) (23) 



and 



Jf~ 1 »(*-(")>-' I (x(0) + g(0)) 3 ( M(0) i m\ (24 ) 

The friction 7(0) and mass parameters M(0) are expressed in terms of first and second 
derivatives of the response function at u — 0. Their modification to the multi dimensional 
case are of the forms 

W oh-^ Lo ^,,_„ (25) 



w=0 



and 



M, 



^(0) = g 



l<3 2 i>(u;) 



&u 2 



w=0 



2 &u 2 



w=0 



(26) 



For obvious reasons, expressions (f2~5|) , ( 26 ) may be referred to as the "zero frequency limit". 



3 The inclusion of pairing degrees of freedom 

In this section we are going to address pair correlations. This shall be done in two steps. 
At first we concentrate on the pairing mode alone to subsequently address the general 
case of pairing plus a shape degree of freedom. 



3.1 The conventional pairing model 

The generator for the pairing mode may be defined as the following pairing field operator 
pt 

^ = £44- (27) 
k 

The a\ and a k are the creation and annihilation operators for the normal vacuum. The 
summation over single-particle states k is meant not to include the time reversed ones, 
which are denoted by k. Commonly one does not start with the associated mean field 
Hamiltonian, as we did in the first section, but with a separable two body interaction 
instead, which may be constructed from the P and P'. The corresponding many body 
Hamiltonian may thus be written as 

H v = H sm - GP ] P (28) 

where the coefficient G of the two body interaction, the coupling constant, stands for a 
constant pairing matrix. The first part H sm is meant to represent motion of independent 
particles, and hence may be written as 

H sm = £ t k {a\a k + a\a- k ) (29) 
k 

with the e k being the single-particle energies. How this H sm relates to the shape degrees 
of freedom shall be of no concern for the moment. As indicated earlier this question will 
be addressed later. 

To establish relations to the discussion in the previous section let us introduce the 
mean field approximation to fl28|) by writing 

H v = H sm - A(pt + p) + |_ + H res (30) 

with 

A = G(P) = G(P j ) (31) 

and 

Hres = -Gj2( a k a l ~ ( a l a l))( a j a j ~ ( a j a j)) ( 32 ) 

kj 



Nothing else has been done but rewriting (]28 ) with the help of the c-number variable A. 



Actually at this stage it may be unclear yet whether or not the averages of P and P^ are 
identical. In the worst case we might simply introduce a complex A. However, as we shall 
see a real one will do. 

Evidently, the mean field part of (^) may be identified as 

H? mf = H sm - A(pt + P) + ^ (33) 

Actually, in spirit of the last section the word "renormalized" is appropriate here; it is 
justified because of the c-number term on the very right. Indeed, this term is necessary in 
order to have the expectation value of H^ m f represent the total energy, see e.g. sect. 3. 1.5 
of ||, in particular eq.(3.1.45). For the present model this c-number term is nothing else 
but the Sp t of (|1|). In summary, the mean field approximation to ( |28|) simply implies to 



neglect the residual interaction H res ; see also @ or |i5|| . 

One might be inclined to associate the H res of (^) with the Vjg} introduced in (g). 
This is not the appropriate connection, however, as the V^} is meant to simulate collisional 
damping. There can be little doubt that it would be asking too much if one were to deduce 
the latter from a separable interaction of the type given in (^). Indeed, later on we are 
going to take into account this latter type of residual interaction, albeit in special form, 
namely through self-energies, which in some sense may be considered a generalized mean 
field approximation. 

The A introduced in ( PH| ) will henceforth be considered the collective variable for the 
pairing degree of freedom and the associated mean field Hamiltonian will be the one of 
(|3"3"|). The alert reader may have noticed some similarity of ([H]) with flT3"|), indicating a 
possible connection of the G to the coupling constant k introduced in the previous section; 
this latter relation will be worked out in more detail below. Notice, however, an essential 
difference between (|13|) and fl3~l"|). Whereas the former case refers to a time dependent 
situation in which the expectation values have a definite meaning. Those appearing in 
(|3"T|), on the contrary, are not yet fully specified. Indeed, the Hamiltonian given in ( |3"BD 
is of practical use only if the A may at first be treated as a free parameter — as it ought 
to be for a collective variable — discarding this subsidiary condition ft3l\). It is only after 
the A has been assigned a special value (which for the static case will be that of global 
equilibrium), or a well defined function of time (as in the dynamic case discussed above 
for the general situation, with its specification to the case of pairing to come below) that 
these expectation values will be defined properly. 

3.1.1 Transformation to quasi-particles 

In the following we are going to diagonalize the Hamiltonian H^mf by transforming to 
quasi-particles in strict sense. In doing so we will not make use of the subsidiary condition 



(31). The relations to be discussed now will hold true in the more general sense, not just 
at equilibrium which we will look at below. Actually it is at this place that we exploit 
the fact of A being a free parameter. 

When accounting for pair correlations it is difficult to conserve particle number. For 
this reason we will refer to the grand canonical ensemble introduced in @ and (H). This 
means to replace the H^ m f by 

8%, = H? mf ~ l*N (34) 



fixing the Lagrange multiplier /i by requiring particle number conservation on average, 
viz 

(N) = N with N = J2( a W + a hk) (35) 



A' 



To diagonalize the Hamiltonian H^f ^ us perform first the Bogolyubov transforma- 
tion 

a k = u k a k + v k a\ a k = u k a k - v k a\ (36) 

to new operators a k , a\ . . . which themselves are required to satisfy the anti-commutation 
rules of Fermions. As one knows, this implies the constraints 

u\ + v 2 k = 1. (37) 

Next we write 



where 



H rmf = U o + H' n + #20 ( 38 ) 

A 2 

k k u 

H 'ii = J2i( e k - ^)( u l - v l) + 2Au k v k )](a\a k + a\a k ) , 

k 

H 2o = J2K e k - ^)2u k v k - A{u\ - vl)](ala\ + a k a k ) (39) 



following common notation. To get from (|3q)-(|39|) the independent quasi-particle Hamil- 
tonian we must chose the coefficients u k and v k of the transformation (|36|) to make the 
last term vanish: if 20 = 0- As ^ * s seen fr° m (§1); this condition has the form 

(e k - fi)2u k v k - A(u 2 k - v 2 k ) = (40) 

This equation, together with normalization condition (|37|), can be solved as 



2 1 / ., , e k~ fi \ 2 1 / 1 £k ~ V I -in 

Mt= H 1+ vfa-^+W ' Vt = H 1 " 75 - t? ~ -* ! ( ' 

Finally, the Hj^f becomes 

H^ f = U' Q + £ E k (ala k + a\a k ) (42) 
k 

with 



E k = J(e k -fiy + Ai. (43) 



3.1.2 Thermodynamics of quasi particles 



Having the quasi-particle Hamiltonian H^f °f (42) at our disposal we may now look at 
the density operator for the quasi-static equilibrium for the pairing degrees of freedom. 
It is a special case of (|) it reads 

p qs (A, T) = i eM-H^f/T) Z> = tr exp(-^' / /T) (44) 



(Please observe that the A still is to be considered a free parameter). For a Hamiltonian 
like the one of (^) representing independent motion of quasi-particles it is not difficult 
to calculate various quantities of interest. For the grand potential one gets 

fi = -TlogZ' = -n-E k )- 2T$>g(l + exp(-E k /T)) + , (45) 

k k u 

and for the particle number 

on ^ 

N = - ir = 2'£n k 46 

where 

2\ T 1 (-, e k - I 1 , i E k 
n k =\u k u k /- = v k -ryu k -i " 

with 



n k = (a{a k r = v\ + (ul - vl)n T k = - ( 1 - tanh , (47) 



= (4«fc) qs = (1 + exp(E fc /T))- 1 (48) 

being the probability that the state is excited. Let V be the average value of the pairing 
field operator for the ensemble flip); one finds 



7> = <P) qs = (P^ 8 = (49) 

k 

with 

fc = (44) qs = u k v k (l - 2nJ) = — (1 - 2r£). (50) 
The internal energy may be written as 

E = (H^ f r + »N = (i^ m/ > qs = 2 £ **n* - 2AP + ^ (51) 

and for the free energy one gets 

A 2 

T = Q. + UN = 2Y,£kn k -2KP + — -TS (52) 
k G 

with the entropy S being given by 

S = 2^s k with s k = [n T k lognl + (l-n T k )\og(l-nl)] (53) 
k 

At various places we need the derivatives with respect to A rather than the energies 
themselves. In principle, they could be obtained from ([51]) and ([52]) . It is much easier, 
however, to exploit eq.(^). For the pairing degree of freedom, and for the present purpose, 
the appropriate Hamiltonian to be used there is the H^ m t of (j33|) . Obviously, its derivative 
is given by 

dHfJA) /A ... 2A 

dA =-( P + pt )+-G ( 54 ) 

Thus (0) leads to 

(5Q |M0) -.(«^) (55) 



It is reassuring to see that the same result may be obtained starting from the expression 
for the free energy given in indeed. 
Obviously, this force vanishes at 

Gr 

thus defining the value A eq at which either the internal energy at fixed entropy or the free 
energy at fixed temperature become minimal. Notice that (|56|) is an implicit equation for 
A cq with the functional form of V(A) being given by together with (|50| ) and the 
expressions given before for Ef. as well as for the thermal occupation number: Eq . (p6|) is 
nothing else but the gap equation at finite temperature, which can easily be brought to 
the conventional form 

2 _^ l-2n£ _^ tanh(E fc /(2T)) (5?) 
G k Ek k Ek 



As one may recognize, the relation fl56|) is of the form given in (pT|), which in this 
context has obtained a specific meaning. For a static situation, the averages appearing 
in PT|) are those of thermodynamic equilibrium, in which case the subsidiary condition 



(|31|) is valid for A replaced by A eq . As this special form is obtained after minimizing the 
total energy, one may understand why sometimes it is referred to as the "self-consistency 
condition" (for the separable interaction we started with). 

3.1.3 Remarks on the coupling constant 

In flU]) a coupling constant k has been introduced which depends on the collective variables 
as well as on excitation. Let us see how this k relates to the G of (|28|). A direct application 
of (0) leads to 

ld 2 H v (A)\ qs 9 
-j^ = ( g% ) + (Xpp(o> = 0) - X T pp) = ^ + (Xpp^ = 0) - xIp) (58) 



For the second equation we have simply made use of (p4|). The response function and 
the isothermal susceptibility appearing there may be traced back to (|TT|) and ( |T5| ) if 
only the is replaced by — (P + P^). Within the independent quasi-particle model their 
calculation is straightforward. (The evaluation of the response in the more general context 
will be given below, in particular in the Appendix). As the susceptibility may be written 
as Xpp — ZcfP/dA it may be obtained by differentiating (fI9D with respect to A keeping 
temperature and particle number fixed. One gets 

X T pp = Xpp(0)-2J:^^ (59) 

(The xpp{uS) is the response function calculated for fixed particle number, for details see 
below). Estimating numerically the difference Xpp ~ Xpp( u — 0) from it is seen to 
be negligibly small as compared to the leading term 2/G. Even at T = 2 MeV it makes 
up less than 10%. Exploiting this fact we may conclude 



- k A w G/2 , to put -k A = G/2 



(60) 



Evidently, for this comparison of &a and G, there is no need to worry about the different 
signs or the appearance of the factor of 2, which are simply due to the fact of using in 
(|28|) a form commonly found in the literature. However, a few words might be in order to 
explain why this relation has turned out not to be an exact one. First of all, we may recall, 
that at T = this would be so in strict sense. Secondly, it had been mentioned already 
when discussing the appearance of the k in (|13|), that for the nucleus as an isolated system 
a more appropriate form would be the one in which in ([14]) the isothermal susceptibility 
gets replaced by the adiabatic one x ad defined for fixed entropy. This difference, x ad— x(0), 
is known to vanish identically for ergodic systems, see ||, in which case —k& equals G/2 
exactly, indeed. Thirdly, it can be said for diabolic motion of independent quasi-particles 
this difference again would vanish identically. This observation follows from the fact for 
such a model diabatic motion is given for fixed occupation numbers in which case the 
do not vary with the collective variables (like the A in the present case). 



3.2 The general mean field for pairing and shape variables 

After this detour into the peculiarities of the pairing model we may now go back to the 
starting point in Sec. |2| to combine the collective variable A for pairing with those for 
the shape. Actually, to keep the numerical effort on a manageable level, we restricted the 
computations to just one shape variable, the elongation parameter Q = Ri 2 /2Rq which 
describes the distance between the left and right centers of mass. To have a dimension- 
less quantity this distance is commonly divided by the diameter of the sphere of identical 
volume, see ||. 

To construct the Hamiltonian H sm (xi,pi, Q^) of ([I]) for the set = Q, A, we simply 



need to choose for the H sm of C|33|) the Hamiltonian H sm (x i} p i} Q) of the deformed shell 
model. The quasi-static density operator for the grand canonical ensemble is given by 
@. The generators of collective motion defined in (j3|), specialized now to = Q,A, 
turn out to be as follows. The one for the pairing mode is given already by (|54"|). Let us 
write it as 

dH rmf (Q,A) ~, 2A 

Fa = ^ } = ~(P + P f ) + (61) 

It explicitly contains the c- number correction 2A/G. The equivalent form for the shape 
degree of freedom reads 



F, 



dH rmf {Q,A) dH sm (Q,A) 



Q 



88 



pot 



dQ dQ Q=Qo 8Q Q=Qo dQ 



= F- 



d£ 



pot 



(62) 



where the last term in the end would have to be calculated through the Strutinsky renor- 
malization, see below. It should be noted however, that these corrections may be omitted 
when calculating the response functions. They simply drop out of the commutators in 

Next we turn to the coupling constants (|14D. Again, following the discussion of ( pCf ) the 
one for the pairing degree is solely determined by G, in that we have — k& = G/2 = —1/k&. 
The one for the shape is given by 

dlF(Q,A,T) 

- «0 = + Xff{uj = 0) (63) 



This one will depend both on our collective variables Q and A as well as on temperature 
T. In principle, there could be a non-diagonal one. Following strictly the formulas given 
above we would find 

d 2 F(Q, A, T) T 
~ KQA = OQdA ^fp\ u = °) = Xfp ~ Xfp{uj = 0) (64) 



The expression on the very right of (64) follows from ( j55|) once its Q-dependence is 
taken into account. We find it wiser, however, to simply put this quantity equal to zero, 
essentially for two reasons, beyond that of simplicity: 1) Like before for AA case, the 
difference of the susceptibilities Xfp ~ Xfp(^ = 0) can be expected to be small at small 
T. 2) Again, we are faced with the fact that the shape and the pairing degrees of freedom 
are not really based on the same footing — which by the way is a feature common to 
almost all microscopic approaches. On the one hand, for the pairing degree of freedom 
there is no mean field at our disposal which might be considered as realistic as that for 
the shape degrees of freedom (completed by the Strutinsky renormalization). As a matter 
of fact, it was constructed from a separable two body interaction. 3) Conversely, for the 
shape degrees of freedom it would not make much sense to start from a similar schematic 
force. We simply believe our approximating the mean field in this case to be better. 

3.2.1 Strutinsky renormalization 

Let us address the calculation of the free energy J 7 , needed for instance for the coupling 
constant k,q given in (|63|) . One might exploit the expression (|52 ) derived within inde- 



pendent quasi-particles model. As is well known, such an approximation holds true only 
in the neighborhood of the Fermi level. But such deficiencies can easily be overcome by 
applying the Strutinsky renormalization procedure |2B|, 24] . Thus we represent T as the 
sum of the liquid drop part and the shell correction 

T = T LDM + (65) 
Unfortunately, we can not use here the standard expression for the shell correction to the 



ground state energy given in |24[. The derivation of this expression rely very much on 
the gap equation which does not hold away from equilibrium with respect to A. For an 
arbitrary A we should proceed in slightly different way. Namely, we may make use of the 
advantage of Strutinsky's energy theorem according to which the shell correction to the 
total energy of the system of interacting particles can be obtained within a mean field 
approximation. In our case it means that we can use the free energy derived not with 



the Hamiltonian ( 28|) but with independent quasi-particles Hamiltonian (f^). Namely, we 



will define the shell correction to the free energy as 

8 J 7 — T — T (66) 

where J 7 is given by fl52|) and T is its average. 

The numerical computation of the average free energy T is done in two steps. Firstly, 
in eqs.(|52|) we substitute summation over the single-particle states by integration over the 
single-particle energies with their density of states being that of the independent particle 
model _ 

/oo 
g s (e)de where g s (e) = ^ 8{e k - e) (67) 
-oo , 



In a second step we replace the density g s (e) by the smoothed quantity 



9s(e) -> g(e) = — / g s (e')f av (- J de' = — J^/™ [— — - 

7av J-oo \ 7av / 7av ^ \ ^av , 



(6* 



where f av (x) is the smoothing function of the shell correction method, see [23, 24]. In 
this way the smoothed free energy becomes 



A 2 

k k k 



(69) 



The averaged quantities n k , s k and (j) k are given by (see also |25 

n k = n(ejfe) = — / n(e)f av [ - — — ) de 

7av J-oo \ 7av / 



~ ~ 1 f°° ( 

s k = s(e k ) = — / s(e)f av 

7av J-oo \ 



e - e k 

7av , 



de 



(e*) = — / 0(e) /«, ( - — — J cte 

7av -/-oo V 7av / 



(70) 



With these quantities at our disposal we may express the first order shell correction bT 
as 

6F{A, T) = 2 $> fe 5n fc - A<50 fe - TSs k ) (71) 
k 

with 5n k = n k — n k , etc. In the vicinity of the Fermi energy the quantities n k , <p k and 
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Figure 1: The shell correction to the free energy ([71]) for neutrons in the spherical Woods- 
Saxon potential as function of temperature for A = (left-hand-side) and as function of 
A for T = (right-hand-side). The neutron numbers are indicated in the figure. The 
parameters of the Strutinsky smoothing procedure are: averaging interval 7 av = 8 MeV, 
order of the curvature polynomial: M — 6. 



s k are strongly varying functions of e k , but they behave rather smoothly far away from 
it. As a consequence, the 8n k , 5<p k and Ss k differ substantially from zero only for states 
close to the Fermi energy This implies that the sums in ( |7lD are rather insensitive to 
a variation of their limits. In particular, the problem of logarithmic convergence in the 
solving gap equation does not appear in (|7"i~[). Thus one need not restrict the summation 
in ([TIP to the so called pairing window. 

Fig.|l| shows the dependence of the shell correction (|7l| ) as function of temperature for 
A = and as function of A for T = 0. It is seen that SJ r (A,T) decreases both with 
temperature as well as with A. The feature may be understood as follows. With growing 
T or A, the occupation numbers n k of the states close to the Fermi energy become more 
and more diffuse and thus resemble more the average value n k . The same holds true for 
the s k and <fi k . Consequently, the absolute values of the Sn k etc. decrease and, hence, the 
sum ([H]) becomes smaller. The temperature at which the shell structure disappears is 
approximately equal to 1.5 MeV. Here, the shell correction makes up only a few percent 
of its value at T = 0. The decrease of the shell correction with the pairing gap (at fixed 
temperature) is less rapid. Even for A = 2 MeV the shell correction still amounts to a 
few dozen percent of its value at A = 0. 



3.2.2 Intrinsic response function for independent quasi-particles 

The time-dependent /z^-response functions for "intrinsic" or nucleonic motion are to be 
calculated starting from the definition (|ll|). To begin with we first address the model of 
independent quasi-particles, discarding for the moment the modifications necessary for 
"collisional damping", which will be discussed below. Moreover, we want to concentrate 
on the FF function. The results for the FP and the PP response are given in the 
Appendix. 

In the quasi-particle representation the (single-particle) operator F can be written as 
F = Y, FkkZvl + F kjikj (otlotj + a\a 3 ) + F kjVkj (4 a j + a 3 a k ) (72) 

k kj kj 

where 

Vkj = u k Vj + v k Uj , S, kj = u k Uj - v k Vj (73) 



Inserting (|72|) into (|TTD after somewhat lengthy but straightforward calculations one gets 



XF F (t) = - nJ))e k3 \F k3 \ 2 sm((E k - E 3 )t/h) 

kj 

29{t) 



J>£ + nj - l) V 2 kj \F kj \ 2 sin((E k + Ej)t/h) (74) 



h kj 

Notice please, that in the first sum no diagonal components of the -terms survive. In 
the end this is due to the fact that the corresponding terms of the operator F commute 
with the Hamiltonian (RSI) and do not contribute to the response function ([73]). That is 
why the first sum in (|74|) is marked by a prime. In the second sum of ([T4]) both diagonal 
and non-diagonal components contribute. The Fourier transform of ([74]) lead to 



(75) 



As it is seen from the last equation, the first line of ( [75] ) does not contribute to the response 
function for T = since in this case all n£ = 0. For nonzero temperature both sums of 
( [75] ) do contribute. Please observe that the second line only represents contributions from 
excitations above 2A, as E k + Ej > 2A. Conversely, as Ej. — Ej may be arbitrarily small, 
contributions to the truly low frequency region may only come from the first line in ([75|). 
Apparently they may built up only at finite temperatures. 
Eq.(|74[) can be brought to the more compact form 

Xff® = E + E £K< " ^Oe^vl^l 2 sm ( E *°- E i*' \ t (76) 

k=j sj^s' kj^j ss' 

by introducing both positive as well as negative quasi-particle energies 



E ks = s^{t k - /i) 2 + A 2 , with 8 = ±1 (77) 
together with the corresponding amplitudes Uk s and v ks 

1 ft e k-H \ 1/2 s ( e fc -At y/ 2 

Uto= ^l 1 --^rJ ' Vh - = ^{ 1 + ^r) (78) 



The quantity £ksjs' is defined similarly to ([73]) 

£,ksjs' = UksUj s ' — VksVjs' (79) 



With the help of the quantities (|77| ) and (79) the Fourier transform of (|76|) can be put 
into a form similar to the one at A = 0, 

xffM = ^-{E E+EE)4V ^ w~f ) + J F *\ 2 ^ 

2h k=j s^s' k^j ss' huJ - K E ks ~ E js .) + IE 

In the same way one may obtain the expressions for the other components of the intrinsic 
response tensor, see the Appendix. 

3.2.3 Intrinsic response functions for collisional damping 

To account for collisional widths of the quasi-particle states we proceed similarly to the 
no paring case |], §j. Let us explain details again at the example of the FF-response 
function. The essence of the procedure is most easily seen after first re-writing ( |80| ) in 
the form 

xf F {u) = -^-{E E + EEKL,vl^l 2 x 

k=j s^=s' kj^j ss' 

dnn T (Q) \g ks (fl)gfJ{Q - oo - is) + g js/ (Q)gi°J (Q + u + is)} (81) 
where Q^J (oj ± is) is the Green function for independent quasi-particles 



The spectral density g ks is related to G k ]{oJ ± is) by 

g ks (aj) = i(G ( kJ(u + is) - g£>(u, - is)) (83) 

As argued in [^] the Green function G k ]{w ± is) must be modified by introducing 
a spreading of the quasi-particle width. This may be achieved by introducing into the 
Green function a complex self-energy E(uj ± is, A, T) = £'(u;, A, T) =)= iT(uj, A, T)/2 

^ (u, ± is) -> £ fcs (u, ± fe) = 1 (84) 

ftu; — L ks — 2j (a;) ± u (co>)/2 

The way of how £ is to be calculated in the presence of pairing will be discussed in the 
next subsection. With Gks an d Q ks given the response function (^1|) takes the form 

Xff(uj) = {EE+E E}tlsjs>Xk sj Au)\F kj \ 2 (85) 

k=j s^fis 1 k^j ss' 

where the amplitudes Xkaja'i^) are defined in the same way as for the system without 
pairing, see M. 

/oo 
——n T (VL) [Q ks (VL)g jsl (VL -u-ie) + Q js >(£l)Gks(tt + u + is)} (86) 
-oo Inn 

The only difference is that the energies t k — or e 3 - — /i in eq.(4.2) of || are replaced by 
Eks or Ejv 

The folding integral in (|86|) is computed by means of the residue theorem closing 
the integration contour in the lower half of the frequency plane. The residues which 
contribute to the integral (jS6|) are that of Gks{Q ± is) (mind (|33"D) and n T (Q). In case of 
the frequency independent width T = T(A, A, T) used here (see also sect. [D| below) the 
poles of QksiQ ± ie) are given by 

Ml ± = E ks =p zT/2 (87) 

and the poles of n T (f2) are given by the Matsubara frequencies 

Ml„ = ±z?rT(2n + 1), n = 0,1,2,... (88) 

The computation of the folding integral (BH) with a frequency dependent width r(u;,T) 
without pairing is discussed in detail in M. 



3.3 The quasi-particle strength function 

The Green function Q ks (K3I) contains the frequency dependent self-energy S. Let us first 
focus our attention on its imaginary part —ImE = T/2, which we will estimate from 
the collision term of the Boltzmann equation for the scattering of two quasi-particles. Its 
standard form can be found in the review of Baym and Pethik |26] for ordinary particles 
and in |27], ZS| for Bogolyubov quasi-particles. As is well known, for A = and for small 
excitations the width (decay rate) Ydihuj) can be written as 



T d (hu J ) = h(huJ-^ 2 + 7r 2 T 2 ] 
1 o 



(89) 



For Bogolyubov quasi-particles no analytic expression for this decay rate Td(fiu) is avail- 
able, rather this quantity has to be computed numerically. In the present work we followed 
the publications ^2| where T^(tiu) was expressed as 



T d (huj) = ^ r J J J de 2 de 3 de 4 5(huj+£2+S 3 +S4)[n^nl+(l-n^)(l-n^)(l-nl} (90) 

Here, the nj represent the thermal equilibrium distribution n T (Si) = (1 + exp(^/T)) _1 
for the quasi-particles, the energies of which are defined as 

S k = E k ^^- (91) 

They are positive (negative) for states arising from single particle states above (below) the 
Fermi energy. The T depends on the averaged matrix element of the residual interaction, 
the effective mass m* at the Fermi energy and a cut-off momentum q c : 

1 / m* \ 3 q c 



T \2nh) /)/.- 



(34|flr|12> (92) 



In deriving ( [90|) an expansion in the momenta around the Fermi momentum was used, 
and the q c is an upper limit for the momentum transfer involved in the scattering, which 



is assumed to be small compared to the Fermi momentum g c (< pf), see [26|. Therefore 



the result (190) is valid for small excitations only. Fortunately, it is this region around the 
Fermi energy which gives the main contribution to the response functions. In case is 
needed also for larger energies a correction becomes necessary. This problem is not specific 
to pairing, it already exists for the zero gap limit where it was solved by introducing one 



more parameter, see [22, 29]. Generalizing this procedure to the case of pair correlations, 



we may assume the final form of decay width to be given by 

r = TT^T? (93) 

with the cut-off parameter c (in energy) chosen to be independent of A and equal to the 
value c = 20 MeV given already in |22[| . Likewise, for the Tq of ( |9~2"D we take the same 
constant as before, namely T = 33 MeV (see p2]j ). 

The computation of the folding integral ([H]) with the frequency dependent width T 
requires the knowledge of T in the whole range of u from zero to infinity. This turns 
out to be a rather difficult task. Already a calculation of the through Q90D ) would 
involve a double integral, which in the end would be too time consuming. To reduce the 
computation time to manageable level we chose to work with a constant T, taken at the 
Fermi energy. This is to say that in the present numerical evaluations we replaced T(hu)) 
by T = T(hou = A). This quantity is shown in Fig.[| as function of temperature for a 
few values of the pairing gap. As might be expected on general grounds, for T = the 



collisional width (^) is zero for all hu < 2A (incidentally, in [22] it was found to vanish 
in the regime huj < 3A). It is also seen that V is the smaller the larger is A. This effect 
is especially strong for small temperatures, say T < 0.5 MeV. 

Once the r(hu, A, T) is known the real part of S is given by Kramers- Kronig relation. 
For systems without pairing it was shown analytically [RJ that E' ~ hu — fi. The same is 
true in the presence of pairing. The numerical results exhibit the £' to vanish at the Fermi 
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Figure 2: The collisional width (^) taken at the Fermi energy (hu = A) as function of the 
temperature. Different curves correspond to different values of pairing gap A indicated 
in the figure. 



energy and to be very small in its vicinity. Within the approximations used here, namely 
to take the width at the Fermi energy where £' = 0, the real part of the self-energy does 
not contribute to the response functions. 

Finally, we should like to mention some shortcomings of the derivation of ( PU| ) given 
in [p2]] , which so far we have not been able to improve, but which for the present works. 
1) In the vicinity of the critical temperature T c , where the pairing gap is small, not 
only the scattering but also the coalescence and decay of quasi-particles should be taken 



into account, see |27j. 2) For small temperatures, where the pairing gap is large, the 
introduction of the energies Sk, instead of quasi-particle energies Ek — which make the 
collision integral resemble that for ordinary particles as closely possible — does not appear 
to be well justified; moreover, one would think that the decay rate should also be influenced 
by the amplitudes u\ and v\. It is very likely, of course, that such finer details will not 
modify our final results too much, as we are using drastic approximations to the width 
anyway. More detailed investigations are required to clarify these points. 



4 Average particle number conservation 

In this chapter we address the problem of particle number conservation. As one knows, 
violation of any symmetry may lead to spurious states. Their existence in the present 
model may be seen with the help of the FF-response function, for instance. Its poles are 
associated with the energies of the systems' excited states. As it is seen from QT5| ) some 
of the poles correspond to pair excited states a]~aj,|0) with energies 2Ek- The number of 
these states is equal to the number of terms in the sums over k or j in (|75|) ), one more 



than the right value, as argued in |L6]. This extra states is the spurious one. Fortunately, 
it is possible to construct response functions where such contributions are removed, as 
shall be explained now. 



4.1 Number conserved response 



4.1.1 Modified intrinsic response 

For the independent quasi-particle model applied to the quasi-static picture, it is the 
chemical potential [i which guarantees the conservation of the average particle number 
determined by eq . fl46|) . However, the dynamical density matrix will differ from the quasi- 
static density p qs . Therefore, in a time dependent picture the average particle number 
(N) t does not automatically coincide with N. Consequently, spurious contributions pro- 
portional to S(N) t might appear. Following [[njand ||I8| , one may correct this drawback 



requiring as a constraint 5(N) t = to hold true, thus warranting that the average particle 
number is conserved also in dynamical sense. 

To implement such a constraint into the theory let us recall first that the Hamiltonian 
(138|) depends on the shape parameter Q, the chemical potential fi and A. The change 
with time of the average of some operator, say (F) t , will generally depend on the variation 
of all three parameters Q, /i and A. To evaluate the deviation of (F) t one may use the 
set of equations (|T2|) which in the present case read: 

6(F) U = - X ff(uj)5Q(u) + X FN(uj)5fi(u) + X fp(uj)5A(u) (94) 

6((P + Pt)) w = -x pf {lo)5Q{lo) + XpnWh{") + Xpp^)5A{uj) (95) 
A similar relation must hold true for 5{N) W 

= - X nf{uj)SQ{uj) + X nn(u)8h(u) + X np{u)5A(u) (96) 

As mentioned previously, the expressions for the FN-, A^P-response functions may be 
found in the Appendix. The condition 6{N) t = or 6{N) U = together with fl96|) allows 
one to express 5fi(u) in terms of 5Q(u) and 5A(u) 

x i \ Xnf{u) , . xnp(u) , , f s 

5fi(w) = --5Q(uj) —8A(uj) (97) 

Xnn{u) Xnn{u) 

By substituting ( |9"7| ) into (Q) and fl9~5| ) the latter two equations turn into 

5(F) U = -Xff(uj)5Q(lu) + x fp (uj)5A(uj) 
5{(P + Pt)) w = -x PF (co)5Q(co) + xpp{u)5A{u) (98) 

In this way we get new response functions Xff(w), Xfp(u) and x PP (ui) which account for 
the conservation of the average particle number during a dynamical process. They may 
be constructed from the former one by the simple relations: 

- / \ / \ Xfn{u)xnf{u) 

Xff{u) = Xff(u) ? — r (99) 

- / a / \ Xfn(u)Xnp(u) . . 
Xfp\w) = Xfp{uj) p-r (100) 

- ( A i \ Xpn(u)xnp(u) , . 

Xpp{UJ) = XppM 7-r (101) 

XnnM 



4.1.2 Modified collective response 



We may now use the information contained in (p9|)- (|101j) , fl6~3f ), and flBOf) to construct 
the collective response tensor ([18]). The expressions for the component of XcoiiC^) can De 
obtained directly from the form given in (0). One gets 

, coiu \ _ «q[(«a + Xpp)xff - XfpXfp] fmo , 

Xffk^) — 7 — w — \ — l iUZ J 

C«0 + XffH^a + Xpp) - XfpXpf 

xssm = ( ,; g "f: Fp - - - (103) 

C«0 + XffH^a + Xpp) - XfpXpf 

and 

coil/ \ _ «a[(kq + Xpp)Xpp - XppXpp] (mA \ 

Xpp{U) — 7 — ; T7 — ; r l iU4 J 

C«Q + XppJI^a + Xpp) - XfpXpf 

The collective Q- vibrational frequency (Qi) is defined by the lowest pole of collective 
response function ( |102|) or by the equation 

(kq + Xff{^i)){k a + Xpp(fii)) - XFp(fii)xpp(^i) = 0. (105) 



It is not difficult to convince oneself that ( |105| ) may as well be interpreted as the secular 
equation of the system fl94T)-(|96|), provided one obeys the subsidiary condition 5{N) U = 
as well as the Q- and A-components of self-consistency conditions (0). In case we are 
interested in vibrations around a minimum of the potential energy, both with respect to 
Q and A, we may identify Q m with Q , A m with A . Thus the Q- and A-components of 
q(uj) will coincide with 8Q(u) or 5A(u) and fljjD will attain the form 

5(F) uj = k q 8Q{uj), *<P + pt) w = |*A(u;) (106) 



Excluding 8(F) u , 5(P + P^) W and 5(iV) w from the system (0)-|§) with the help of flTUBD 
together with 8(N) U = one will get the following system of equations 

- (xff{m) + k q )5Q(uj) +xfn{u)8/j,(uj) +xfp{u)8A(uj) = 

-Xnf(u})5Q(u) +xnn(u)8ii(lu) +xnp(u)8A(u) = (107) 
- X pf(oo)8Q(uj) +xpn{oj)8^{u) +(xpp(^) ~ 2/G)5A(lu) = 

The secular equations is then given by putting the determinant of the matrix behind this 
set equal to zero. Indeed, this leads to the same result as in ( |105| ). Below, this fact shall 
be exploited further for special cases. 



4.2 RPA at zero thermal excitations 



It should be worth while to compare our results with those obtained in previous treatments 
at zero temperature. Using the expressions for the intrinsic response functions given in 
the Appendix for T = one derives from (|107| ) the following secular equation 
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(108) 



This equation may be compared with the secular equations obtained in [19|, |20[ and used 
for investigations of vibrational states in medium and heavy nuclei. The determinants 
of [|l9j and [|20| differ from each other by a constant multiplier so that the collective 
frequencies are the same. The one obtained in [B(J (let us denote it as DetRPA(cu)) writes 
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(109) 

To establish a relation between ( |108j ) and ( p.09| ) let us first recall that in [ 19], pU[ vibrations 
around the minimum of the potential energy were considered. In this case, because of 
the linearization involved, the pairing gap may be defined through the solution of the 
gap equation (corresponding to the minimum) and one may replace 1/G by ^ k \j2E k . 
Consequently, one has 
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(110) 



Considering this relation, the determinants ( |108|) and ( |109|) look very similar. Indeed, by 
explicit calculation one may verify the following relation 
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It is clear from ( |1 1 1| ) that our secular equation ( |108| ) defines the nontrivial part of the 
secular equation of [19|, pD[ . The spurious mode corresponding to u = is not present 
among the solutions of (|108|) , which are the same as those found in flTj^, f20|j . 

Later on we are going to evaluate transport coefficients. They are largely determined 
by the strength distributions associated to the frequency dependence of response functions, 
like xff{u) and x ff {uj). Our numerical computations show that the shift in the strength 
distribution due to correction terms like Xfn(^)xfn(^) / Xnn(^) is small. This implies 
that particle number conservation has only little influence on the transport coefficients. 
The correction terms mostly cause shifts in the positions of the peaks of (the imaginary 
parts of) the response function, whereas the widths of the latter remain almost unchanged. 
As a consequence, the friction coefficients deduced from xff(u) and x ff (uj) are almost 
identical. Not much difference is seen is also in the mass parameter and the local stiffness. 



4.2.1 Pairing without shape oscillations 

For a comparison with the standard pairing theory let us consider the particular case 
where the coupling between Q-mode and A-mode may be neglected (x F p(u) = 0). In 
this case the secular equation ( |105| ) separates into individual ones for the Q- and the 
A-mode. The latter can be written as 



(112) 



Replacing 2/G in the same way as above by J2k^/^k and using the explicit expressions 
for the response functions, eq. (|112|) is transformed to 



<^-- 2 »(?W^V)) -(?5§Ffe) (U3) 

After taking the square root eq. (|113|) is reduced to 



x _ J% 2 u 2 -4A 2 ±2(e k -a) , s 

V ww^A -° (114) 

This equation is identical to that obtained in [in| within quasi-boson approximation. 
In this sense essential results of the theory of vibrational states in paired systems are 
reproduced by our linear response approach. 



5 Numerical results 

In this section we are going to present numerical results for response functions and trans- 
port coefficients. We concentrate on deformations along the fission path of the nucleus 
22A Th, which has been under investigation before. Indeed, present one is an extension of 
previous publications where temperatures above T = lMeV had been considered 



with pairing discarded. Like in || we use the independent particle Hamiltonian based on 
the deformed Woods-Saxon potential with the nuclear shape parameterized in terms of 



Cassini ovaloids [31]. The Cassini ovaloids are defined by rotating the curve 



p(z, e) = #0 



a 4 + Aez 2 /Rl - z 2 /R 2 - e 



1/2 



(115) 



around the z-axis, with z and p being cylindrical coordinates. The constant a is deter- 
mined from volume conservation, implying that the family of shapes ( |115|) depends only 
on one deformation parameter e. As is easily recognized from ( |115|) the value of e = 
corresponds to a sphere. For < e < 0.4 the form ( |115|) resembles very much that of a 
spheroid with the ratio of the axes given by 

shorter axes 1— 2e/3 

= — (116) 

longer axes 1 + e/3 

At e ~ 0.5 a neck appears and at e = 1.0 the nucleus separates into two fragments. Like 
in ||, instead of e we will use the parameter Q = Ri 2 /2Rq, which measures the distance 
R12 between the left and right center of masses divided by the diameter 2Rq of the sphere 
of identical volume. Besides having a simple physical meaning, the Q = R\2/2Rq allows 
one to relate to other shape parameterizations of deformed potentials. 

In principle, the pairing gap is different for neutrons and protons. Thus, for quan- 
titative studies one should introduce gaps for protons and neutrons as two independent 
dynamical variables. However, in the present work we would like to concentrate more 
on questions of principle nature, as a general investigating of the effects of pairing on 
transport coefficients. Thus we feel justified to compute the response functions putting 
A p = A™ = A and to use this A as a collective parameter common for protons and 
neutrons. This will simplify also the presentation of our results. 
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Figure 3: The frequency dependence of the imaginary part of the intrinsic response func- 
tions (|^)-( |T0Tl ). The computations were done with the spherical Woods-Saxon potential 
for the nucleus 224 Th. The values of the temperature T and A are indicated in the fig- 
ure. The dashed and solid curves correspond to different choices of the collisional widths, 
namely ( [39|) and (|90|) , respectively (together with (l93|)). 



5.1 Information from the intrinsic response 

We begin looking at the intrinsic response functions x^u, as the most basic quantity from 
which all local transport properties may be deduced, once the static energy is known. 
For instance, it allows for an immediate calculation of inertia and friction within the zero 
frequency limit, which for not too small temperature may give reasonably accurate results, 
in particular in the version suggested in ||. Moreover, this response function is a basic 
element in the construction of the collective response function, which in turn allows one 
to evaluate the implication of self-consistency on transport coefficients. It is thus of great 
interest to look first at the effects pairing has on the Xw m Fig-ll we show the dissipative 
parts of the susceptibilities given in (|99|)-( |10lD as function of frequency, for T = 0.5MeV 
and two values of A, for spherical configurations. The functions presented in all three 
parts of the figure exhibit a narrow peak structure. These peaks represent the excitations 
the system would have if there were no correlations due to collective motion. 

Let us first look at the effects the collisional width T(hu = A, A, T) has on width and 
height of the peaks in the response functions. To clarify this influence, and consequently 
the one on the transport coefficients, in Fig.|3| two sets of curves are plotted: One computed 
with T(hu = A, A = 0,T) (dashed curves) and another one with the T(hu = A, A, T) 
given by eqs.fl93|) and ( j9Dp . As it is seen from the figure, the peaks computed with a finite 
A are higher and more narrow as compared with those for the case A = 0. This is to 
be expected since, as it is seen from Fig.[|for fixed temperature, T(hu = A, A, T) is the 
smaller the larger is A. For fixed A the width T(A, A, T) increases with temperature and, 
as the computations demonstrate, the response functions become smoother. 

Next we examine the influence of the spectrum of the individual excitations. As 
compared to the unpaired case, the response functions shown in Fig.^| exhibit a totally 
different structure: The lowest peak of considerable strength occurs at or above the 
frequency hu « 2A, and the strength in the region tkv < 2A is very small. In case of the 
PP- response, for instance, the position of the peaks is defined by the zeros of eq. 
The two peaks seen at the bottom part of Fiq.^| represent the contributions from the lowest 
solutions of eq. ( |114j ) for neutrons and protons. For a spherically symmetric potential the 
lowest peaks are well isolated from the rest. This is because the positions of the peaks 
of the PP-response function is mainly determined by the energy 2Ek of pair excitations. 
This distance is especially large for spherical configurations, simply because of the large 
degeneracies of the single-particle states. For deformed shapes the distance to the next 
peaks becomes much smaller. 

Obviously, these features will have great impact on transport coefficients. This is 
most apparent for the friction coefficients defined in zero frequency limit: They will be 
the smaller the larger is A. This property also shows up more or less strongly for the self- 
consistent calculation. Let us therefore analyze the structure of the response functions at 
small excitations in more detail. It is easily recognized that it may be traced back to the 
various contributions to formula ( ff5|) (or the corresponding one resulting after collisions 
are taken into account). As we will see, the contributions from the two lines behave very 
differently, depending very sensitively on temperature. Let first look at the case of finite 
temperatures. Here, the first line of (f75|), resulting from the terms with — Ej\ < 2 A, 
will have finite contributions even in the region hu < 2A. However, on the scale used 
in the figure, they are very small. Fig.^ shows the (dissipative part of the) response 
function in the interval below 2A in more detail, exhibiting clearly the existence of such 



contributions. 
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Figure 4: The frequency dependence of the imaginary parts of the response (solid curves) 
and correlation (dotted curves) functions. The computations are done with the spherical 
(left) and deformed (right, e = 0.2) potentials. 



Since around u = 0, the region which is of special interest, this dissipative part has to 
vanish exactly, it is more instructive to use the correlation function tj)''(u}), related to the 
former by the fluctuation dissipation theorem (for the nucleonic degrees of freedom), see 
0. One has to replace the commutator of (|11|) by an anti-commutator and must account 
for the unperturbed averages relative to which the deviations 5F^ = F M — of the 
operators are to be defined, 



^(t) = \(([H(t) } F v ] + )- w<^» 



(117) 



Calculating the averages in ( |117|) in the same way as for the response functions, for 
independent quasi-particles the Fourier transform of the FF-correlation function turns 
out to be 
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with E k j = Ek ± Ej. The account of collisional damping is carried out in full analogy 
to the response function. The ip' FF {uj) is shown in Fig.[| by the dotted line. Comparing 
with the response function (full line) the properties just mentioned are clearly visible. 
In addition, another important and generic feature is revealed. At certain deformations 
the correlation function not only is finite at and around uj = but shows a more or 
less pronounced and sharp peak. For the special case shown here, which corresponds to 
a deformation of e = 0.2, it appears to lie exactly at to = 0. In this sense it reminds 
one of the "heat pole", which in previous publications was seen to become more relevant 



at larger temperatures, see [|32| and ||. However, there are cases where such a peak 
is shifted to small but finite frequencies. These peaks appear whenever a pair of levels 
close to the Fermi energy approach each other; in the present case their distance is about 
\e k — 6j\ ~ 0.02 MeV which is definitely smaller than the width associated with such 
excitations. 

Having the correlation function at one's disposal, the friction coefficient in zero fre- 
quency limit can be expressed as 

7ff(0) = ^ f (uj = 0)/(2T) (119) 

Hence, for this limit it becomes apparent how peaks in the correlation function at u = 
are related to peaks in the deformation dependence of the friction coefficient. For the 
special example discussed here it may be said that the above mentioned peak in ip FF at 
oj = results in a peak of the friction coefficient shown in Fig.|5], namely at Ri2/2R « 0.43 
which corresponds to e = 0.2. 
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Figure 5: The zero frequency limit (p5|) for the friction coefficient in the shape degree of 
freedom Q as function of Q for T = 0.5 MeV and A = 1.0 MeV. The right bottom and 
top left parts of figures show separately the contribution to the friction of nondiagonal 
Fkj and diagonal matrix elements. The bottom right part contains the stiffness of 
the static energy at T = 0.5 MeV. The curves with stars mark the value obtained by 
averaging over deformation on the interval A(Ri2/2Rq) = 0.08. 

Let us turn to the extreme case of T = for which A will be finite and, more important, 
for which the collisional width will be zero at least for hu up to Hu = 2A; recall the 
discussion below eq. (p3|) . Look at the forms ([75]) and ( |118|) the response and correlation 



functions attain for the case of no collisions. Then, only contributions from the excitations 
above 2A survive. In case of the presence of collisions, in the region uj < 2A both the 
response and correlation functions may differ from zero only if tails from distant peaks 
extend into this regime. This is very unlikely to happen if the convolution integral in 
(|86|) is calculated exactly with the proper frequency dependent widths. Hence, friction 
in zero frequency limit will be strictly equal to zero. However, in practical computations 
where the convolution integrals in fl86|) are computed with constant T's such tails might 
be present, which in turn may lead to a finite though small value of the friction coefficient. 
The choice used in our calculations, namely T = T(A, A, T), as described below (p6[), does 
not suffer from this drawback. 

As the zero frequency limit provides the simplest definition of friction and inertia (see 
(25) and (^)) respectively, we like to stick to this version for a moment. In Fig.[5| we 
exhibit the deformation dependence of friction for the Q-mode, as evaluated from (P5|). 
In this figure the rapid oscillations with deformation we spoke of can be seen. The inertia 
looks similar to friction. The formulas of the zero frequency limit (|25| ) - (|26|) are simple 
enough to allow one separating the contributions from diagonal and non-diagonal parts of 
the response function (|[5|) (diagonal here means the contribution with k = j). The source 
of fluctuations of these two parts turns out to be different. The diagonal sum (contribution 
of the diagonal part of the second sum in (|75|)) is affected strongly by the density of states 
near the Fermi energy, which is typical of the appearance of shell effects, of such a type 
as they also appear in the static energy. They exist as long as the temperature does 
not exceed a typical value of the order of 1 — 2 MeV. The stiffness C (0) of the static 
energy is shown in the right bottom part of Fig.[|. Indeed, the same kind of fluctuations 
appear, only anti-correlated to the diagonal component of friction. The fluctuations of the 
non-diagonal part are caused by the quasi-crossings of levels near the Fermi energy like 
those discussed in connection to Fig.f|. For each sharp peak in the left-hand-side of Fig.|| 
one can find in the single-particle spectrum one or a few quasi-crossings which contribute 
90% or so to friction or inertia. The height of such peaks is the larger the smaller the 
minimal distance is between quasi-crossing levels. From the analysis of the single-particle 
spectrum computed with the deformed Woods-Saxon potential it follows that some of 
these minimal distances are extremely small - of the order of 0.02 — 0.03MeV. 

For several reasons we claim that to large extent such closely spaced levels are short- 
comings of the underlying shell model, and that they may lead to unphysical consequences. 
First of all, it should be noted that they are intimately related to the choice of the shape 
parameter, in particular to the restriction to only a few of them, more precisely to exactly 
one in the present case. In a truly multi-dimensional landscape such quasi-crossings would 
happen much less frequently as there is always enough room for avoiding them. Secondly, 
residual interactions would reduce their significance further. Indeed, from the statistics 
of experimentally measured nuclear states it is known that levels with small spacings are 
very seldom; the distribution of nearest neighbor spacings is of Wigner type rather than of 
Poisson type. Such interactions would repel the approaching states and in this way would 
reduce the matrix elements with the generator F of collective motion. Evidently, it should 
be the aim of any microscopic theory to deal with each of such problems individually. At 
present, we are not yet able to do so. Rather, we suggest to overcome these difficulties 
in another way, namely by averaging the transport coefficients over deformation. The 
averaging interval should be large enough to smooth out the rapid oscillations, as well 
as the finer details of the shell structure, but at the same time small enough to preserve 
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Figure 6: The averaged values of the components of the friction (left-hand-side) and 
inertial (right-hand-side) tensors obtained within zero frequency limit d25|) ,(|26| ) as function 
of deformation for temperature T = 0.5 MeV. The dotted, dashed and solid curves 
correspond to A = 0.5, 1.0 and 1.5 MeV. 



gross shell effects. This procedure may perhaps be motivated further by recalling the 
following facts, (i) Evidently, any structure in the transport coefficients of greater detail 
cannot be seen experimentally. Usually, one is happy if one succeeds to identify gross shell 
behavior, (ii) Finally, these transport coefficients are to be used in a transport equation of 
Fokker-Planck or Langevin type which account for genuine fluctuations in the collective 
variable. Without any doubt, such fluctuations will smooth out the detailed structure 
in the coefficients in most natural way. Incidentally, we may phrase our problem in a 
different way, namely in attributing it to a deficiency of the mean field approximation, at 
least partly. 

The averaged friction coefficient and stiffness are shown by the curves with stars in 
Fig.[5| As it is seen the frequent oscillations are gone while the typical gross structure 
remains. This latter feature is clearly visible at the case of the stiffness which now shows 
the typical behavior expected for a potential with a second (and perhaps third) minimum. 

In Fig.^| we exhibit the averaged tensors of friction and inertia are shown, for T = 
0.5 MeV and for a few values of the pairing gap A. Both friction and inertia still oscillate 
as function of the deformation, but in much weaker fashion than those shown in Fig.|5]. 
Now they correlate with the fluctuations of the stiffness or the shell correction. The 
non-diagonal terms, namely QA-friction and inertia, oscillate around an average value 
which is very small. Very likely this non-diagonal components of friction or inertia may 
simply be neglected. The AA-friction or inertia fluctuate much less as compared with 
the corresponding quantities in Q. The reason is that the matrix elements of F (which 
are of strongly peaked structure) do not contribute to AA-friction or inertia, see fllOip . 
The quantities which matter here most are the single-particle energies themselves which 
are relatively smooth functions of deformation. As seen from Fig.^], an increase of A has 
two effects, both on friction as well as on inertia: They become smaller and their Q- 
dependence becomes smoother with growing A. In case of the QQ-transport coefficients 
the substantial contribution comes from the diagonal component of the F operator (diag- 
onal term in eq. (|74]) ). For not too small A his diagonal contribution can be estimated to 
behave like oc 1/A 2 both for friction and inertia. The temperature dependence (for fixed 
A) of friction and inertia is the same as in the no pairing case: Whereas friction increases 
with temperature, inertia decreases. 

5.2 Information from the collective response 

In Fig.|7| a few examples of the tensor of the collective response are shown for some values 
of A and temperature. A common feature of all figures is as follows: The collective 
frequency is shifted to smaller frequencies as compared to the lowest peak of the intrinsic 
response. The magnitude of the shift depends on the value of the static response 
relative to the stiffness 0^(0), see Fig. 6 of [|l2j for a demonstration. The smaller is the 
stiffness the larger is the shift of collective strength to lower frequencies. In case of the Q- 
mode the static stiffness is much smaller than x'ff(®) anc ^' hence, the collective strength 
concentrates at low frequencies. The situation turns out very different for the A-mode, 
whose stiffness Caa approximately equals the static response x'pp(O) implying that the 
shift of the collective peak is small. 

In ideal cases the X^O^) * s dominated by one peak in the low frequency region. This 
peak may then be fitted rather accurately by the response function of the damped oscil- 
lator, in which way the self-consistent transport coefficients are defined, see [[|. Unfortu- 
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Figure 7: The frequency dependence of the imaginary part of collective FF- and PP- 
response functions, associated to the intrinsic functions shown in Fig|| The solid and 
dashed curves correspond to (|102|) , (|104|) and to the fit to (|102|) , (|104|) by the response 
function of a damped oscillator. 



nately, in reality one often encounters a different situation in that the collective strength 
splits over several closely placed peaks which may even overlap. This situation is typical 
for a multi-dimensional case, for which it has been demonstrated in JO]] of how transport 
coefficients may still be deduced through a generalized fitting procedure. Indeed, most 
likely the present situation ought also be treated in this way. Unfortunately, the method 
suggested in [O] is tedious and very time consuming. In order to circumvent problems of 
this type we prefer to define transport coefficients through the approximation defined in 
(|2"T|), which still exploits the information contained in the form (kx co11 ~ 1 (uj)K,) flL/ of the 
collective response tensor. The fit to the polynomial of second order given on the right 
hand side of (|2"T|) can easily be done for any function, even if x coll 



-i, 



does not look 

much like a polynomial. In such cases the transport coefficients may be considered as 
some frequency average over the interval of the fit. Unfortunately, averaging in frequency 
does not change much the fluctuations with respect to the deformation. That is why 
both quantities were also averaged over the deformation like before. The results of such 
a procedure are shown in Fig.^| for the diagonal components of friction and inertia. As it 
is seen from the figure, the transport coefficients obtained in this way exhibit fluctuations 
with deformation similar, but noticeably different to those of the zero frequency limit, 
which have been shown before in Fig.|6|. Moreover, the mean value of both components of 
friction and inertia tensor now are larger than before. The difference is especially large in 
case of AA-friction, the reason being that at small u the imaginary part of the response 
function for the pairing degree of freedom xpp{&) increases with uo much slower as com- 
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Figure 8: The components of the friction (left-hand-side) and inertia (right-hand-side) 
tensor obtained by a fit to ( |102| ),( |lO4] ) through the response function of a damped oscillator 
according to (^Tj) as function of the deformation for temperature T = 0.5 MeV. The 
dotted, dashed and solid curves correspond to A = 0.5, 1.0 and 1.5 MeV. 



pared with that of a damped oscillator. Consequently, 7aa(0) is smaller than 7aa- This 
is a clear hint that in such a case the zero frequency limit may not be considered a good 
approximation. 

5.3 Temperature dependence of transport coefficients 

In the past various models have been presented to justify nuclear dissipation, even more 
exist to evaluate the friction coefficient numerically. In [33[] a compilation of data has 



been put together in which theory is confronted with experimental evidence. The latter 
is obtained by comparing solutions of the "macroscopic equations" (of Fokker-Planck or 
Langevin type) with experimentally observable quantities, the input into such equations 
being chosen phenomenologically. It was seen that the predictions which the theoretical 
models give for the effective, local collective width T = 7/M of the fission mode - 
sometimes referred to as the "reduced friction coefficient" (3 — deviate by as much as two 
orders of magnitude. It can be said that the results of |53J classified there by "linear 



response theory" (LRT) were in good agreement with those required by experimental 



evidence. As compared to the application of LRT mentioned in [Q, new calculations 
have in the mean time been presented in || ffOfl . Grossly speaking their results lie in the 
same regime as those of |34j]. 

As all such microscopic calculations depend on some quantities which are not known 



all too precisely, it has been argued previously (see e.g. [§, gOf) that one should look for the 



temperature dependence of the effective transport coefficients. Indeed, here the various 
theories differ from one another in very pronounced ways. The wall formula, for instance, 
predicts a friction coefficient which is practically independent of T. As demonstrated in 
p2| , p, such a picture is obtained as the macroscopic limit of LRT, provided the latter 
is applied to purely independent particle motion, viz if all influences of collisions are 
discarded. Truth is that the latter have important implications on the T-dependence 
of friction. Two body viscosity, for instance, comes about in the regime of "collisional 
dominance", in which viscosity is known to decrease with temperature, in ideal cases 
like 1/T 2 . Although the notion "viscosity" intrinsically refers to a nuclear liquid, in the 
case of genuine " collisional dominance" such a T-dependence also prevails for the friction 
coefficient of the low frequency collective modes of finite nuclei, see e.g. [[JIJ. A similar 
dependence was found in and p7| . It is true that these models do not rely on a 



hydrodynamic description, but in each one collisions are assumed to dominate collective 
dynamics in one way or another. In |3^] and it has been discussed how this behavior 
of damping may be obtained within LRT. 

Fortunately, a possible T-dependence of nuclear dissipation has also been studied by 



interpreting the experimental results, see [pq, 39, HQ . Here we shall not pursue any 



details of these studies. It may suffice to summarize the general claims put forward in 
these publications, which is that dissipation increases with T at small excitations, leveling 
off or even decreasing at larger temperatures. Indeed, generally speaking, such behaviour 
was found in our previous calculations [ 32| , at least qualitatively. How large friction may 
become at intermediate temperature and how quickly it decreases afterwards is still an 
unsettled question. It depends both on the contributions from the heat pole as well as on 
the approximations in which collisions are treated, for details see |52 and ||. As stressed 
earlier, in the present paper we mainly address the regime of very small temperatures 
where pairing must not be neglected. Unfortunately, to the best of our knowledge, in this 
regime, neither experimental information is available, nor to date is there any theoretical 
model with which our results might be compared. Our predictions are as follows. 

The ones for the fission mode are displayed in Fig.[5| To simplify matters the pairing 
mode has been treated in equilibrium, or "local" equilibrium, rather, if the latter is being 
defined with respect to the collective variable Q. In other words, at given shape the gap 
equation (|57|) has been solved for various temperatures. The resulting A's then have been 
used for the calculation of the transport coefficients for the Q-mode displayed in Fig|| 
For the system considered pair correlations disappear above T ~ 0.5 MeV. Indeed, it is 
only below such temperatures that the transport coefficients deviate from the unpaired 
case. At first sight the deviation may look small. A closer look, however, shows that 
important modifications of the nature of the transport process are present. Neglecting any 
quantum corrections (for collective dynamics) one realizes that for smaller temperatures 
nuclear dissipation may get so weak that one reaches Kramers' so called "low viscosity" 
limit". This is very interesting on recalling that in nuclear physics experimental results 
are commonly interpreted on the basis of Kramers' simple formula for the "high viscosity 
limit" . Moreover, it so happens that in this regime quantum effects in collective dynamics 
are present. They lead to additional modifications of the decay rate. Unfortunately, at 
present one is able to account for these effects only above some critical temperature, which 
turns out to be less than the T where pairing disappears. There is no room to elaborate 



further on details of collective dynamics. Instead we like to refer to a recent Letter [|T 
as well as to forthcoming papers. 




T (MeV) T (MeV) T (MeV) 



Figure 9: The average QQ-friction coefficient 7qq, reduced friction coefficient (3qq = 
^qq/Mqq and the damping factor r] QQ = ^ Q q/2^CqqM Q q at the ground state (top) and 

fission barrier (bottom) of 224 Th as function of temperature. The dashed curves show the 
results obtained neglecting pairing. 



Let us turn to the pairing mode A now, for which in Fig. 10 we show the damping 
factor t?aa = 7aa/2-\/Caa^aa as function of temperature for various values of A. Here, 
the latter is to be understood as a collective variable, not fixed by the gap equation. The 
values presented in this figure represent averages along the fission valley. The //aa is seen 
to be very small, implying the corresponding pairing vibrations to be strongly under- 
damped. Most likely this means that on the way from saddle to scission the pairing mode 
does not equilibrize. The consequences such an observation may have on odd-even effects 
in mass distributions etc. will have to be the subject of some further studies. 



6 Summary 

A detailed study of transport coefficients has presented for temperatures below T = lMeV 
which finally aim at a description of the dynamics of fissioning nuclei. Pairing effects 
were accounted for within the independent quasi-particle approximation. An appropriate 
treatment of particle number conservation has been achieved by modifying the response 
functions of the locally harmonic approximation. For harmonic vibrations without damp- 
ing it has been demonstrated that the same secular equation for the collective frequencies 
is obtained as those derived earlier within common RPA. In our version, the pairing gap A 
is introduced as an independent dynamical variable, similar to the parameters specifying 
the shape of nuclear surface. The response functions and the tensors of collective friction 
and inertia have been calculated for a realistic, deformed Woods-Saxon potential. The 
components of friction and inertia tensors were examined as function of the pairing gap 
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Figure 10: The damping parameter of pairing vibrations tjaa = 7aa/2v / Caa^aa as the 
function of temperature for few fixed values of A. The damping parameter is averaged 
along the fission valley of nucleus 224 Th. 



and the deformation of the nucleus. With respect to the latter strong variations were 
found of all transport coefficients. They are caused both by shell effects as well as by 
avoided crossings of single-particle levels. It has been argued that such variations may 
be considered as carrying unphysical features, for which reason they have been smoothed 
out by averaging over the shape in such a way as to keep gross shell effects. As could be 
expected, the friction coefficient decreases with increasing A. Finally we may say that the 
temperature dependence of the friction coefficient obtained in our model is in qualitative 
agreement with the conclusions reached in |I8|, [39|, |4(| , namely that dissipation increases 
with T at small excitations. 
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A Intrinsic response functions 

The time dependent response function XFFif) is defined by ([TT|) . To obtain XPNit) one 
should replace P 7 (t), F I (s) in (|Tlf) by (P + P^) 1 (t) , N I (s). In the same way one can obtain 
expressions for the other response functions, too. For the operators P and N one may 
use the quasiparticle representation like (|72|) 

P + P ] = 2u k v k ~ ^u k v k {a\a k + a\a- k ) + Y,( u t ~ v l)( a W k + a ~k a k) 

k k k 

N = 2v l + J2( u l - v k)(°i a k + a\a k ) + 2u k v k (ala\ + a k a k ) 

k k k 

(120) 



The final expressions for FP or FiV-response functions turn out to be much simpler than 
the one presented in fl5U| ) since they contain only diagonal sums (k = j). 



X fn(oj) = ]T(2r^-i)(2 Mfc t; fc ) 2 F fc 



1 



1 



Xfp(uj) = J2( 2n k ~ l)2u k v k {u 2 k - v 2 k )F, 

k 

XNN(uj) = J2( 2n k ~ l)(2«fcVfc) 

k 

XNP{UJ) = J2( 2n k - l)2Mfcffc(Mfc - v 2 k ) 



hu — 2E k + ie hu + 2E k + ie 
1 1 



kk 



hco — 2E k + ie hto + 2E k + 
1 1 



hu — 2E k + ie hu + 2E k + ie 
1 1 



X pp(oj) = EC 2 "* -!)(«*- 



,2\2 



ftu; — 2F& + ^u; + 2Ffc + ie 
1 1 



hu — 2E k + ie hu + 2E k + ie 



(121) 



The account for the collisional damping for FP- FN- and other response function is 
carried out in the same way as for FF-response and leads to 



Xfn(u 
Xfp{u 
Xnn(uj 

Xnp{u 
Xpp(uj 



Y^{^u k v k ) 2 F kk E Xksks>{u) 

k sj^s' 

2u k v k (u 2 k - v 2 k )F kk E Xksks'(uj) 

k SytzS' 

J2(2u k V k ) 2 E Xksks'i^) 

2u k V k (u 2 k - V 2 k ) E Xksks'(uj) 
k s^s' 

Y,( u l- V l) 2 E Xksksi^j) 

k sj^s' 



(122) 



with Xksks'{u) given by 
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